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Abstract 

We systematically study the dissociation of H2 molecules on the /3-Ga203 (100)B surface, with 
the influences of surface oxygen vacancy being considered. After introducing the surface oxygen 
vacancy, the nearest topmost 0(1) atom becomes very active, and hydrogen molecules become 
much easier to dissociate. 
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I. INTRODUCTION 



In ambient conditions, the stable form of gallium oxide is the monoclinic structure of 
/3-Ga 2 3 which belongs to the group of transparent oxides with a large band gap of 

4.8 eV, and exhibits both conduction and luminescence properties [4H3]. The /3-Ga 2 3 is 
generally an n-type semiconductor due to oxygen vacancies individually compensated by 
two electrons forming shallow donors js-10]. However, the number of oxygen vacancies 
present in the material depends on the growth atmosphere, and so the electrical character 



11- 



Because 



of the compound can vary, in a tunable way, from insulating to conductive 
of its tunable optical, magnetic and electrical properties, the fl-G^O^ has attracted much 
research interest for many_technological applications jlO ]. For example, Ga 2 3 is useful in 



the fabrication of masors field-effect devices I I], switching memory devices 
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mii, and 



transmitting ultraviolet light [18|, |19|. By monitoring its electronic properties at different 



chemical atmosphere, the /3-Ga 2 03 can also be used as sensors to detect hydrogen gases 



20- 



23] . To detect hydrogen molecules, it is revealed that oxygen vacancies are needed on the 



Ga 2 03 surface to explain the stretching frequencies of Ga-H bonds observed in the infrared 
spectroscopy experiment 24j. However, the special mechanism for oxygen vacancies to help 



dissociate hydrogen molecules has not been discussed yet. 

In addition, the H 2 -Ga 2 03 interaction is also important to some catalytic reactions. Cat- 
alysts containing supported gallium are known to be active in light alkane dehydrogenation 



and aromatization 



24 



27] . Recently, it is highlighted that /3-Ga 2 3 alone is able to dissoci- 



ate hydrogen molecules at temperatures higher than 500 K and then hydrogenate adsorbed 



28|. Although 



29 



;he oxidation state of gallium 



30] , the role of the geometric 



C0 2 stepwise from formate to methoxy groups 
ions has been of a matter of debate in catalytic reactions 
environment or surface coordinations has not been emphasized. 

There are two types of coordination for Ga ions in the /3-Ga 2 03, namely tetrahedral and 
octahedral, and referred to as Ga(I) and Ga(II). There are also three types of oxygen ions, 
referred to as 0(1), O(II) and O(III). Oxygen 0(1) and O(II) lie respectively in and out of 
the symmetry plane (see Fig. 1(a)). They are both in threefold coordination while O(III) 



is in fourfold coordination |6 
(100) plane 



23 
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10fl . The most frequent cleavage plane for /3-Ga 2 3 is the 



31] . A recent theoretical study of the (100) surface presents two possible 



surface terminations, A and B 



23, 



31] . Termination A is characterized by rows of oxygen, 
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and termination B by rows of nearest neighbor gallium and oxygen. In our study, we mainly 
focus on the surface relaxations for the stoichiometric and oxygen vacancy included /3-Ga2C>3 
(100)B surfaces, and the dissociation prorperties of hydrogen molecules on them. The rest 
of the paper is organized as follows. In Sec. II, we describe our first-principles calculation 
methods; in Sec. Ill, we present in detail our calculated results, including the comparisons 
of the electronic structures of the /3-Ga203 (100)B surface including oxygen vacancy or not, 
and the dissociation of hydrogen molecules on the surfaces; Finally, the conclusions are given 
in Sec. IV. 



II. COMPUTATIONAL METHOD 



Our calculations are performed within density functional theory (DFT) using the Vi 



enna ab-initio simulation package (VASP) |32j. The PW91 [33j generalized gradient ap- 
proximation and the projector-augmented wave potential js^J are employed to describe the 
exchange-correlation energy and the electron-ion interaction, respectively. The cutoff en- 
ergy for the plane wave expansion is set to 400 eV, which is large enough to make the error 
from calculations of the adsorption energy below 0.01 eV. The stoichiometric and oxygen 
vacancy included /3-Ga2C>3 (100)B surfaces are modeled by a slab composed of 30 atomic 
layers and a vacuum region of 20 A. The 1x2 supercell with the size of 5.93 A x 6.11 A is 
used to study the H 2 adsorption. Our test calculations have shown that this supercell is 
sufficiently large to avoid the interaction between adjacent hydrogen molecules. Integration 



over the Brillouin zone is done using the Monkhorst-Pack scheme [35( with 7x7x1 point 
grids. A Fermi broadening 36[ of 0.05 eV is chosen to smear the occupation of the bands 
around the Fermi energy (Ej) by a finite-T Fermi function and extrapolating to T = K. 
During geometry optimizations, the bottom 10 layers of the /3-Ga203 (100)B surface are 
fixed while other Ga and O atoms are free to relax until the forces on them are less than 
0.01 eV/A. The calculation of the potential energy surface for molecular H2 is interpolated 
to 350 points with different bond length (g?h-h) and molecular height (/iH 2 ) of H 2 at each 
surface site. The calculated lattice constants of the /3-Ga203 crystal are a=12.45 A, 6=3.05 
A, and c=5.93 A, and #=103.9°, in good agreement with the experimental and theoretically 



calculated values 
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23]. The calculated bond length of the free H 2 molecule is 0.75 A, also 



in good agreement with the experimental value of 0.74 A [371 ] . 
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III. RESULTS AND DISCUSSION 



First, we do geometry optimizations for the stoichiometric /3-Ga203 (100) surface. There 
are two different types of (100) surface terminations, respectively called as the ft-G^O^ 



(100)A and (100)B surfaces 23j, |31[. The (100)A surface is formed by cleaving the Ga- 



O(II) bonds above the O(II) atoms while leaving the below Ga-O(II) bonds intact, and 
terminated in rows of 0(II)'s lying along the [010] direction. The (100)B surface is formed 
by cleaving the Ga-O(II) bonds below the O(II) atoms, and terminated in nearest-neighbor 
rows of Ga(II) and O(III) atoms, each singly unsaturated. After geometry optimizations, 
we find that the topmost 0(1) atoms move out on the fl-G^O?, (100)B surface, and form 
0(I)-Ga(II)-0(II) rows as shown in Fig. 1(a). The calculated surface energies for the relaxed 
f3-Ga.20s (100)A and (100)B surfaces are respectively 1.19 and 0.76 J-m -2 , according well 



with other ab initio results of 1.13 and 0.68 J-m -2 23] . The surface energy values indicate 
that the fl-G^Oz (100)B surface is much more stable than the f^-G^Os (100)A surface. 

There are two kinds of oxygen vacancies on the relaxed /3-Ga203 (100)B surface, re- 
spectively the 0(1) and O(III) vacancies. After geometry optimizations, we find that the 
/3-Ga 2 03 (100)B surface with the O(III) vacancy is 0.56 eV lower in total energy than with 
the 0(1) vacancy. The surface geometry of the /3-Ga203 (100)B surface with an O(III) va- 
cancy (Ov-surface) is shown in Fig. 1(c), where the possible adsorption sites near the O(III) 
vacancy are depicted. On the Ov-surface, the 0(1) atom has a relatively large distortion. It 
moves 0.28 A along the [001] direction during surface relaxation. 

Besides of introducing geometric distortions, the oxygen vacancy also changes the surface 
electronic structures. Figure 2 shows the surface charge density of the stoichiometric (5- 
Ga203 (100)B surface and the Ov-surface. The contour map shows the charge density 
distribution in the plane that is parallel to and 1.50 A away from the /3-Ga 2 03 (100)B 
surface. And the values are in the unit of e/A 3 . It is clearly shown that the 0(1) atom 
nearest to the O(III) vacancy has more electron distribution than on the stoichiometric 0- 



Ga 2 3 ( 
analysis 



00)B surface, and so does the nearest O(III) atoms. Using the Bader topological 
38| for the charge density, it is found that the nearest 0(1) atom gets 0.28 more 
electrons on the Ov-surface than on the stoichiometric /3-Ga 2 03 (100)B surface, it is also 
found that the two nearest Ga(II) atoms and the nearest Ga(I) atom respectively get 0.61, 
0.62 and 0.34 more electrons on the Ov-surface. These extra electrons are because that the 
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vacancy of a surface O(III) atom supplies 2 more electrons to its nearby atoms. And the 
other 0.15 electrons is around the other surface oxygen atoms. The extra electrons around 
the surface Ga atoms are easy to form Ga-H bonds when hydrogen atoms are introduced. 

The projected density of states (PDOS) for the nearest surface 0(1) atom on the Ov- 
surface is shown in Fig. 3, together with the PDOS for a surface 0(1) atom on the stoichio- 
metric /3-Ga203 (100)B surface. From detailed wavefunction analysis, we find that there are 
strong s-p hybridizations in the electronic states for Ga(I) and Ga(II) atoms on the /3-Ga20 3 
(100)B surface, but the electronic states of surface 0(1) and O(III) atoms do not hybridize 
at all. In the PDOS of an 0(1) atom on the stoichiometric /3-Ga 2 3 (100)B surface, the 
s states lie around 15 eV below the p states, and the p electronic states contribute the 
majority part of the electronic states around the Fermi energy. As shown in Fig. 3, in the 
PDOS of the nearest 0(1) atom on the Ov-surface, there is a new peak (peak 1) below the 
Fermi energy, which represents for the electronic states of the 0.28 extra electrons due to 
the O(III) vacancy. Since electrons near the Fermi energy are easier to join in electronic 
hybridizations, the extra electrons around the surface 0(1) atom make it more active than 
on the stoichiometric Ga203 surface. Moreover, it is found from wavefunction analysis that 
the electronic states of peak 1 tend to occupy the oxygen p z orbital. These facts indicate 
that the electronic states around the 0(1) atoms are very easy to interact with a parallel 
hydrogen molecule. As we will see later, hydrogen molecules adsorbing around the surface 
0(1) atom on the Ov-surface do dissociate more easily. 

After analyzing the surface electronic structures, we then study the adsorption and dis- 
sociation of hydrogen molecules, both on the stoichiometric /3-Ga 2 3 (100)B surface and 
the Ov-surface. As shown in Fig. 1(b), there are five high symmetry sites on the relaxed 
/9-Ga203 (100)B surface, respectively on top of the Ga(I), Ga(II) and oxygen atoms, and in 
the middle of two Ga(I) atoms and two oxygen atoms. In following discussions, we call them 
as the TGI, TG2, TO, BG and BO sites. At each surface site, there are three possible ori- 
entations for an adsorbing hydrogen molecule, respectively along the x ([001]), y ([010]) and 
z ([100]) directions. After employing these notations, we will use TGl-x,t/,z, TG2-x,y, z, 
TO-x, y, z, BG-x, y, z, and BO-x, y, z to represent the considered adsorption channels for 
hydrogen molecules. While studying the adsorption and dissociation of hydrogen molecules 
on the Ov-surface, we choose TO to be the site directly on top of the O(III) vacancy, and 
TGI and TG2 sites to be the sites on top of the nearest Ga(I) and Ga(II) atoms. 



5 



To study the adsorption properties of H 2 molecules, we do geometry optimizations after 
initially putting an H 2 molecule 4.00 A away from the stoichiometric /3-Ga 2 3 (100)B 
surface and the Ov-surface along all the above mentioned high symmetry channels. During 
the geometry optimizations, the bottom 10 layers of the /3-Ga 2 3 (100)B surface are fixed 
while other H, Ga and O atoms are free to relax until the forces on them are less than 0.01 
eV/A. From these relaxation calculations, we find that there are no chemisorption states 
for H 2 molecules on both the stoichiometric /3-Ga 2 3 (100)B surface and the Ov-surface. 
All the H 2 molecules will finally evolve into the physisorption states, with the adsorption 
energies of several tens of meV. 

We calculate the two-dimensional (2D) potential energy surface (PES) cuts to evaluate 
the dissociation properties of hydrogen molecules. Along each adsorption channel, hn 2 ranges 
between 0.40 A and 4.00 A, and o?h-h ranges between 0.45 A and 2.40 A. From the 
calculated 2D PES cuts for H 2 on the stoichiometric /3-Ga 2 3 (100)B surface, we see no 
direct dissociation. As examples, the 2D PES cuts along the TGl-x, y, TG2-x and TO-y 
channels are listed in Figs. 4(a)-(d). As shown in Figs. 4(a)-(c), on the stoichiometric 
/3-Ga 2 3 (100)B surface, hydrogen molecules are very hard to get close to surface Ga atoms. 
The energy barrier for a hydrogen molecule in its molecular length of 0.75 A to be 1.00 
A from the surface Ga(I) and Ga(II) atoms is over than 12 eV. However, hydrogen molecules 
might be easier to get close to the O atoms on the /3-Ga 2 3 (100)B surface. The energy 
barrier for a hydrogen molecule in its molecular length to be 1.00 A from the surface 0(1) 
atom is 5.2 eV, as shown in Fig. 4(d). In total, hydrogen molecules can not dissociate on the 
stoichiometric /3-Ga 2 3 (100)B surface. Even the energy barriers for hydrogen to get closer 
to the surface is huge. These results indicate the inactivity of the stoichiometric /3-Ga 2 3 
(100)B surface to hydrogen molecules. 

Similar to what we do for the 2D PES cuts of hydrogen molecules on the stoichiometric 
/3-Ga 2 3 (100)B surface, we also calculate the 2D PES cuts of hydrogen molecules on the 
Ov-surface. To make a comparison, we list the 2D PES cuts along the TGl-x, ?/, TG2-x 
and TO-y channels in Figs. 5(a)-(d). One can see from Figs. 5(a)-(c) that the adsorption 
of hydrogen molecules does not change too much on top of Ga(I) and Ga(II) atoms on the 
Ov-surface, comparing with that on the stoichiometric /3-Ga 2 3 (100)B surface. However, 
the adsorption of hydrogen molecules change a lot on top of the 0(1) atom. As shown in Fig. 
5(d), there is a clear dissociative adsorption state for hydrogen on the Ov-surface, in which 
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the bond length and molecular height of H 2 are respectively 1.80 and 0.50 A. The dissociation 
energy barrier for a hydrogen molecule to reach this adsorption state is 2.97 eV, which is 
much smaller than the approaching energy barrier on the stoichiometric /3-Ga 2 3 (100)B 
surface, indicating a big enhancement of the surface reactivity to hydrogen molecules. The 
PES results also approves our previous hypothesis that due to the extra electrons from the 
surface O(III) vacancy, the nearby 0(1) atom are easier to interact with hydrogen molecules. 
We also notice that the dissociation energy barrier of hydrogen molecules is still very large. 
But more importantly, we point out that surface oxygen vacancies can greatly enhance 
the activity of the /3-Ga 2 3 (100)B surface to interact with hydrogen. Given more surface 
oxygen vacancies, and an environment with larger hydrogen enthalpy, hydrogen molecules 
can be further easier to dissociate on the /3-Ga 2 3 surface. 

IV. CONCLUSION 

In summary, we have systematically studied the electronic structures of the stoichiometric 
and oxygen vacancy included /3-Ga 2 3 (100)B surfaces, and the dissociation properties of 
hydrogen molecules on them. It is found that due to the oxygen vacancy, the neighboring 
surface Ga and O atoms get some extra electrons than on the stoichiometric Ga 2 3 surface. 
The extra electrons around the surface 0(1) atom are very near to the Fermi energy, and 
distribute on the p z orbital, which are thus very active to hydrogen molecules. By calculating 
the 2D PES cuts for hydrogen molecules on the stoichiometric and oxygen vacancy included 
/9-Ga 2 3 (100)B surfaces, we find that hydrogen molecules are very hard to get close to 
the stoichiometric Ga 2 3 surface, but much easier to dissociate on the oxygen vacancy 
included Ga 2 3 surface. Our results indicate the enhancement of surface activity to hydrogen 
molecules by introducing surface oxygen vacancies. 
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List of captions 



Fig.l (color online) The /3-Ga 2 3 (100)B surface viewed along the (a) [010] and (b) 
[100] directions, and (c) the /3-Ga 2 3 (100)B surface with an O(III) vacancy viewed along 
the [100] direction. Grey and red balls respectively represent Ga and O atoms. The black 
squares in (b) and (c) represent the surface unit cells in our calculations. 

Fig. 2 (Color online) Contour map of the charge density distribution on the (a) stoichio- 
metric and (b) oxygen vacancy included /3-Ga 2 3 (100)B surfaces. The mapping plane is 
1.5 A from the Ga 2 3 (100)B surface. The back squares in dashed lines represent the 
surface unit cells. 

Fig. 3 (Color online) The projected density of states for the 0(1) atom on the stoichio- 
metric and oxygen vacancy included /3-Ga 2 3 (100)B surfaces. Fermi energy is set to be 
the energy zero. 

Fig. 4 The 2D PES cuts for the adsorption of hydrogen molecules along the (a) TGl-x, (b) 
TGI-?/, (c) TG2-a; and (d) TO-y channels on the stoichiometric /3-Ga 2 3 (100)B surface. 
The total energy of a free H 2 molecule plus that of the /3-Ga 2 3 surface is set to be the 
energy zero. 

Fig. 5 The 2D PES cuts for the adsorption of hydrogen molecules along the (a) TGl-x, 
(b) TGl-y, (c) TG2-x and (d) TO-y channels on the oxygen vacancy included /3-Ga 2 3 
(100)B surface. The total energy of a free H 2 molecule plus that of the oxygen vacancy 
included /3-Ga 2 3 surface is set to be the energy zero. 
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FIG. 2: 



12 



O(I) near the O(III) vacancy 
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